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Abstract. Bose-Einstein condensates are described in a mean- field approach by the 
nonlinear Gross-Pitaevskii equation and exhibit phenomena of nonlinear dynamics. 
The eigenstates can undergo bifurcations in such a way that two or more eigenvalues 
and the corresponding wave functions coalesce at critical values of external parameters. 
E.g. in condensates without long-range interactions a stable and an unstable state are 
created in a tangent bifurcation when the scattering length of the contact interaction 
is varied. At the critical point the coalescing states show the properties of an 
exceptional point. In dipolar condensates fingerprints of a pitchfork bifurcation have 
been discovered by Rau et al [Phys. Rev. A, 81:031605(R), 2010]. We present a 
method to uncover all states participating in a pitchfork bifurcation, and investigate 
in detail the signatures of exceptional points related to bifurcations in dipolar 
condensates. For the perturbation by two parameters, viz. the scattering length and a 
parameter breaking the cylindrical symmetry of the harmonic trap, two cases leading 
to different characteristic eigenvalue and eigenvector patterns under cyclic variations 
of the parameters need to be distinguished. The observed structures resemble those 
of three coalescing eigenfunctions obtained by Demange and Graefe [J. Phys. A, 
45:025303, 2012] using perturbation theory for non-Hermitian operators in a linear 
model. Furthermore, the splitting of the exceptional point under symmetry breaking 
in either two or three branching singularities is examined. Characteristic features are 
observed when one, two, or three exceptional points are encircled simultaneously. 



PACS numbers: 03.75.Kk, 03.65.Vf 
1. Introduction 

The macroscopic occupation of the bosonic ground state of ultracold quantum gases 
has been predicted by Bose and Einstein, and at least since the first realisation of 
Bose-Einstein condensates (BECs) [l}]3) they are a central part of experimental and 
theoretical atomic physics. The condensates are typically held in a harmonic trap. The 
particles in the trap interact via short-range contact interactions determined by the s- 
wave scattering length a, which can be varied experimentally with the help of Feshbach 
resonances. In addition long-range interactions exist, e.g., in dipolar condensates [4j or 
in condensates with an attractive 1/r interaction [sIIg]. 
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In a mean-field approach the effective one-particle wave function of the condensate 
is described by the nonlinear Gross-Pitaevskii equation (GPE). The nonlinearity of the 
GPE allows for a variety of phenomena which are impossible in linear quantum systems 
with Hermitian operators. For example, the extended Gross-Pitaevskii equation for a 
BEG without long-range interactions or with an attractive 1/r interaction has a second 
solution which emerges together with the ground state in a tangent bifurcation [6]. At 
the bifurcation point both states coalesce, i.e., the energies and the wave functions 
are identical. It has been shown |7) that the bifurcation point has the properties 
of an "exceptional point" [8 12 . Such points can appear in systems described by 



non-Hermitian matrices which depend on a two-dimensional parameter space. The 
exceptional points are critical points in the parameter space where both the eigenvalues 
and the eigenvectors of the two states pass through a branch point singularity and 
become identical. There is only one linearly independent eigenvector of the two states 
at an exceptional point. 

Bifurcation scenarios become even more substantial in dipolar condensates. With 
a simple variational approach a tangent bifurcation between the ground state and an 
unstable excited state was observed in [13j , and it was shown that the bifurcation point 
has the properties of an exceptional point. However, using an extended variational 
approach with coupled Gaussians for the condensate wave function the ground state 
is created unstable and only becomes stable at an increased value of the scattering 



length 14 . Although no bifurcating states have been directly computed in 14 , the 



stability change indicates the existence of a pitchfork bifurcation, which means that 
three states coalesce at the bifurcation point. 

The detailed investigation of pitchfork bifurcations and signatures of three 
coalescing eigenfunctions in dipolar condensates is the objective of this paper. We 
present a method which exploits the symmetry breaking of the external harmonic trap 
and allows us to reveal all three states involved in the pitchfork bifurcation. The energies 
and the eigenfunctions of the three states coalesce at the bifurcation point, however, the 
degeneracies are lifted when the scattering length is varied or the axial symmetry of the 
harmonic trap is broken. 

The signatures of three coalescing eigenfunctions have been investigated by 
Demange and Graefe 15 using perturbation theory for a model with linear but non- 
Hermitian operators. It was shown that two types of parameter perturbations need 
to be distinguished, which are related to a cubic root and a square root branching 
singularity at the exceptional point and lead to characteristic patterns under cyclic 
variation of the parameters. With an analytic continuation of the nonlinear GPE in the 
two parameters, viz. the scattering length and a parameter breaking the axial symmetry 
of the harmonic trap, we show that the signatures of the three states of a dipolar BEG 
which coalesce in the pitchfork bifurcation exactly resemble those obtained in the linear 
model 15 with complex non-Hermitian matrices. Both, the cubic root and the square 
root branching singularity are observed when encircling the exceptional point either in 
the symmetry breaking parameter or in the scattering length. Furthermore, we will 
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discuss the behaviour of the eigenvalues and eigenvectors when both control parameters 
are simultaneously varied. 

2. Variational approach to dipolar condensates 

An extended variational ansatz with coupled Gaussian functions has been used by Rau 



et al 16 ,17 to compute numerically accurate solutions of the GPE for condensates with 
dipolar interactions. Contrary to the numerical imaginary time evolution of states on 
grids the variational approach allows for the calculation of not only the ground state but 
also excited states, which is a prerequisite for the investigation of bifurcation scenarios 
and exceptional points. In this section we briefly recapitulate the variational method 
and introduce the analytic continuation of the GPE, which then allows for the encircling 
of exceptional points in the complex parameter space. 

We use the time-dependent GPE in dimensionless and particle number scaled units 
(see (T3|), 

/f^(r, t) = [-A + V]ijir, t) = i^^(r, t) , (1) 
to describe dipolar BECs in the mean-field approximation. The total potential 

V = Vt + V, + Va (2) 
is composed of three parts. The external trap is described by the harmonic potential 

K = 7^ + ly + , (3) 

where the '~fx,y,z determine the strength of the trap in the x, y, and z directions. The 
dipolar condensates are typically held in an axisymmetric trap with 73, = 7y. In order 
to describe the trap frequencies in a way that represents the trap geometry and to 
allow us to change the trap symmetry by modifying only one parameter we will use the 
parametrisation 

7x = 7(1 + ^)^/'A-^/^ ly = l{l + s)-'/'X-'/\ l. = l\''' (4) 

where 7 is the average trap frequency and A the ratio between the trap frequency in z- 
direction and the trap frequencies in the x, |/-plane. The asymmetry parameter s breaks, 
for s 7^ 0, the axial symmetry of the trap, and is one of the two parameters used below 
to examine the signatures of exceptional points. The potential 

K = 87ra|^(r,t)|2 

describes the short-range contact interactions resulting from s-wave scattering between 
particles with a the (scaled) scattering length. The scattering length a is the second 
parameter used to examine the signatures of exceptional points. Finally, the potential 

K. = /dV ^-=^;°y W.(r'.i)P (6) 

describes the dipole-dipole interaction between particles with a magnetic moment, with 
d the angle between the vector r — r' and the 2;-axis along which all dipoles are aligned 
by an external magnetic field. 
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In order to solve the GPE Q we use the time-dependent variational principle 
(TDVP) |18|. The idea of the TDVP is to use an ansatz for the wave function depending 
on the set of variational parameters z and then choose these parameters in such a way 



that the quantity / 
functions 



\\H%p — i-?/'p is minimised. We use a superposition of Gaussian 



N 



N 



k=l k=l 

as an ansatz for the wave function, where the are the complex width parameters of 
the fcth Gaussian in a direction, and the real and imaginary parts of the describe the 
phase and the amplitude of the fcth Gaussian, respectively. The variational parameters 
are combined in the complex vector 

^ = (7^•••,7^^i,•••,^:,^^•••,<,^L•••,^f)ec^^. (s) 

The application of the TDVP (for details of the derivations see p6j) leads to the following 
equations of motion. The time derivatives of the variational parameters are given by 
z = h{z), or in components 

2i (A^ + AI + A^) + 



At 



If, 



with 



' 



a 



x,y,z 



(9) 
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The quantities 

system of equations 

Mv = r . 
The matrix M has the form 



M 



with the submatrices 

(O),, = {gW) 

for k,l = 1, . . . ,N and the operators O - 
components of the right-hand side vector r 
in equation ( 10 ) read 

N 

r'f = J2{9'\fV\g') 



and in equation tej) can be written vector V 



), which is obtained by solving the linear 



(10) 
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(12) 
The 



N 1 
1 ' X 1 yi ' 



N 1 
I ' y ) ' 2 5 



(13) 



fe=i 

for f = l,x,y, z and I = 



1, 



A^. All components of the matrix M and the integrals for 



the vector r in equation (13) are listed in Appendix A The components of the matrix 
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M can be expressed analytically. The vector r is more complicated due to the dipole 
interaction. The components 

{9^\fV^y) with/ = l,x,y,z (14) 

lead to expressions which contain elliptic integrals of the form 

3 /"°° dt 
RDix,y,z) = - (15) 
2 Jo ^/{x + t){y + t){z + ty 



These integrals can be calculated numerically using the Carlson algorithm 19 



The stationary states are obtained as fixed points of the equations of motion ([9]), 
i.e. — i/i = and A'^ = with fi the chemical potential, a = x,y, z, and k = 1, . . . , N. 
The set of the variational parameters 7 can be reduced to 7^^ = 7^= — 7I for k = 2, . . . , N 
by using the facts that the wave function is normalised to = 1 and a global phase 

of the wave function is arbitrary. The conditions for a stationary state then simplify to 

^'=7'-y = 0; k = 2,...,N, 

i^ = ij = i^ = 0; k = l,...,N. (16) 

The root search for the fixed points is performed using the Newton-Raphson algorithm. 
The mean-field energy of a stationary state is given by 

E^f ={^\-A + Vt + |(K + , (17) 

and can be used to examine some properties of the state. However, the typical 
bifurcation behaviour is revealed even more clearly by, e.g. the expectation values w 
of the operator — y"^, viz. 

w={tlj\x^-y'\tlj) , (18) 

and therefore this value is used below in most cases to reveal the bifurcations and to 
analyse the signatures of the exceptional points. For the presentation of the splittings 
we use the distance from their mean values. 



^ i=l 

1 

AWj = - X] ' ^'^^^ 



1=1 

i.e. the mean values of the quantities of all A^s states participating in the bifurcation 
(with Ng = 2 for a tangent bifurcation and Ng = 3 for a pitchfork bifurcation) are 
subtracted. The index j indicates the jth state of the bifurcation. With the new 
values the level-splittings at the bifurcations (which are small compared to the absolute 
values) can be clearly observed. 

The linear stability of stationary states is analysed by calculating the eigenvalues 
of the Jacobi matrix 
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which is obtained with the variational parameters spht into their real and imaginary 
parts, viz. 

z = (Rezi,Imzi, . . . ,Rez4Ar,Imz47v) G M^^ . (22) 

If all eigenvalues are purely imaginary, the state is stable, otherwise it is unstable. The 
eigenvalues occur in pairs ±A. If one pair of eigenvalues has a nonzero real part the 
state is unstable in the directions of the associated eigenvectors in the parameter space. 
If two pairs of eigenvalues have a nonzero real component the state is unstable in four 
directions, and so on. 



2.1. Analytic continuation of the GPE 

The parameter s in equation (|4]), which breaks the axial symmetry of the harmonic 
trap, and the scattering length a in equation ([s]) are the two, physically real, parameters 
used to control the system. To investigate a bifurcation point in equation ^ for the 
signatures of an exceptional point it is necessary to encircle the critical value in the 
complex plane. Therefore, the GPE ([T]) must be continued analytically with respect to 
the parameters s and a in such a way that observables, such as the mean-field energy 
i?mf or the value of can become complex. As the variational parameters z G C^^ are 
complex parameters it is necessary to split the components into their real and imaginary 
parts, and to start with the real parameters z G M®^ introduced in equation (22). The 
real components of the vector z can now be continued analytically as 

Zl = + kz/ (23) 

with the real parameters z^ and zf for Z = 1, . . . , 8A^, and the imaginary unit k defined 
by k^ = — 1. It is important to distinguish k from the imaginary unit i for the following 
reason. The complex conjugate of the wave function '?/'(r, z) in equation ([T]) is obtained 
by the replacement i — )■ — i and thus z — t- z*, however, the sign of k in the complex 
continued vector z in equation (23) must not be changed. Formally, the variational 
parameters z can be written as bicomplex numbers 



Zl = zj^ + i Zl = (z/,1 + k zi^2) + i (^«,3 + k zi 



4 



= zi^i + k zi^2 + i zi^3 - j zi^4 (24) 

with i^ = k^ = —1, = 1, ij = ji = k, jk = kj = i, and ik = ki = — j. The complex 
conjugate wave function in this description has the form 

N 

^*{r) = ^e"(^-^''+^^*J''+^'*^'+^'*) , (25) 
1=1 

withal; = (v4^^ + kA|^2)-i(^a,3 + k^a,4) and 7'* = (7{ + k7^) -i(7^ + k7^). For complex 
parameters s and a the analytically continued equations of motion ^ can now be set 
up for the variational parameters extended to bicomplex numbers. The advantage of 
using bicomplex numbers is that the rules for calculations can be easily implemented in 
computer algorithms, e.g. by introducing a new data type and overloading operators. 
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In particular, Carlson's algorithm 19 for the computation of the elliptic integrals 



REi{x,y, z) in equation (15) can be used with nearly no modifications for real, complex. 



and even bicomplex numbers x, y, and z. The linear set of equations (10) can be solved 
either directly with bicomplex numbers or the components of the bicomplex matrix M, 
and the vectors v and r can be split into two complex numbers to set up an ordinary 
complex linear system of equations of twice the dimension. 



The stationary states are obtained as roots of equations (16) extended to bicomplex 
variational parameters. Note that with the analytic continuation observables such as 
the mean-field energy E^f or the value of w are not necessarily real but can become 
complex valued, e.g. w = Wi + kw2 with the i and j components = = 0. In the 
following we use the notation Kew = wi and Im w = W2 in such cases. 

3. Bifurcations of the stationary states 

We study the bifurcation scenarios of the stationary states of a dipolar BEC in a 
harmonic trap for the case that the average trap frequency 7, the ratio A, and the 
symmetry breaking parameter s in equation ^ are kept constant and the scattering 
length a is varied. A sketch of the bifurcation scenario is given in figure [l](a). The ground 
state and an excited state emerge in a tangent bifurcation (marked Tl in figure [T]^a)) 



at a critical scattering length a = oti 13 , 17 . When the scattering length is increased 



both the ground and the excited state can undergo further bifurcations indicated PI, 
P2, and P3 in figure [T|a). 

3.1. Bifurcations of the excited state 

We start with the investigation of the stationary states using a single Gaussian function 
as a variational ansatz for the wave function. The trap parameters are set to 71 = 34000, 
Ai = 6, and s = for a harmonic potential with axial symmetry, which is in the region 
where dipolar condensates with chromium atoms have been realised experimentally 
(4,20. The bifurcation scenario drafted in figure |l](a) (without P2) emerges. At the 



scattering length oti = —0.019 a stable ground state and an unstable excited state are 
created in a tangent bifurcation (Tl). Both wave functions are axisymmetric around 
the z axis. At the scattering length api = —0.0079 the excited state undergoes a 
pitchfork bifurcation PI, where its stability properties change from unstable for one 
pair of stability eigenvalues to unstable for two pairs of stability eigenvalues. If the 
scattering length is further increased all three states which emerge in the pitchfork 
bifurcation are recombined in an inverse pitchfork bifurcation at ap3 = 0.129. At the 
bifurcation points not only the mean-field energies and the expectation values of the 
operator — but also the wave functions coalesce. 

In figure [2]^a)-(b) the detailed results for the pitchfork bifurcation PI around the 
critical scattering length api = —0.0079 are presented. Shown are both the real 
branches and the complex branches obtained by analytic continuation of the GPE. 
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a) 



b) 




a 



a 



Figure 1. (a) Sketch of the bifurcation scenario of a dipolar EEC in an axisymmetric 
trap. Two states emerge in a tangent bifurcation at Tl. The two branches undergo 
pitchfork bifurcations at PI and P2. The three states created in the pitchfork 
bifurcation PI merge again in an inverse pitchfork bifurcation at P3. (b) Scheme how to 
find the branching states of the pitchfork bifurcation: Starting at s = on the central 
state marked by a cross in the left-hand upper part of the figure the axial symmetry 
of the trap is broken (cross in the lower part of the figure), then the state is followed 
adiabatically by increasing the scattering length, and finally the axial symmetry of the 
trap is restored to end up with a state on one of the symmetry breaking branches in 
the upper part of the figure. 



The characteristic shape of a pitchfork bifurcation is evident in the real part of Aw in 
figure [2]^a), where two new real branches emerge at the bifurcation point api = —0.0079. 
The new emerging states are unstable with one pair of stability eigenvalues with nonzero 
real part, and the corresponding wave functions break the axial symmetry of the trap. 
The two states only differ by a 90° rotation around the z axis. With the analytic 
continuation two additional complex branches are found below the bifurcation point as 
can be seen in the imaginary part of Aw in figure [2]^a). 

The existence of a new branch at a > api is also visible in figure [2]^a) for the splitting 
Aii^mf of the mean-field energy. As the two new states only differ by a 90° rotation around 
the z axis in the axisymmetric trap the mean-field energy of the new branch is twofold 
degenerate. Furthermore, the splitting AE'mf increases quadratically with a — api in 
contrast to the typical square root behaviour of a pitchfork bifurcation, and the mean- 
field energy is real also for those states obtained with the analytic continuation of the 
GPE in the region a < api, i.e. Im (AE'mf) = for all states in figure |2](a). 

When the scattering length is kept constant at a = a-pi and the asymmetry 
parameter s is varied the splittings Ai^mf and Aw show a quite different cubic-root-like 
behaviour as illustrated in figure ^b). The splitting is Ai^^f ~ ■s^/'^ for the mean-field 
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Figure 2. (a)-(b) Level splittings of the mean-field energy iJ^f and the values 
w = {ip\x'^ — y^lV') around the pitchfork bifurcation PI of the excited state. The 
symbols and solid lines mark the results of the GPE using a variational ansatz with a 
single Gaussian function and the linear model introduced in section [3. 3[ respectively. 
In (a) the scattering length is varied for a symmetric trap (s = 0) and in (b) the 
asymmetry parameter s is varied while the scattering length is kept constant (a = api). 
(c)-(d) Results for the pitchfork bifurcation P2 of the ground state. Here, the GPE 
was solved using an ansatz with six coupled Gaussians. 
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energy and Aty ~ s^/^ for the expectation values of the operator — y^, which imphes 
that the sphttings are real-valued only for the central state marked by red plus symbols 
in figure p[b). A more detailed discussion of the splittings based on a linear model with 



non-Hermitian matrices will be presented in section 3.3 



3.2. Bifurcations of the ground state 

Using the simple variational approach with a single Gaussian function {N = 1) the stable 



ground state emerges in a tangent bifurcation 13 . However, the situation becomes more 
complicated when using an improved and extended ansatz with coupled Gaussians for 
the condensate wave function ( I4|[I6[[T7| . For = 6 Gaussian functions and an external 
trap with parameters 72 = 6887, A2 = 7, and s = the ground state is created unstable 
in the tangent bifurcation and then becomes stable at a larger value of the scattering 
length. The stability change exhibits fingerprints of a pitchfork bifurcation, however, the 
numerical search for the new branches emerging in the bifurcation is a very nontrivial 
task in the dimensionally increased parameter space for coupled Gaussians and thus 



have not been computed in 17 . Here we present a method based on breaking the axial 
symmetry of the harmonic trap to obtain the additional states emerging in the pitchfork 
bifurcation. 

The method consists of three steps as illustrated in figure [l]^b). Starting with 
an axisymmetric state (s = 0) at a scattering length below the bifurcation point the 
asymmetry parameter s is varied in a first step to break the axial symmetry of that state 
(s 7^ 0). In a second step the scattering length a is increased to a value greater than 
the critical value of the pitchfork bifurcation. Using small increments of the scattering 
length we can follow the state adiabatically as shown in figure [l](b). In a third step 
the parameter s is changed back to s = for the axisymmetric trap, however, the state 
at the end of this path belongs to one of the new branches emerging in the pitchfork 
bifurcation. The other branch can be reached in the same way with the asymmetry 
parameter s changed in the opposite direction. 

The states of the new branches break the axial symmetry around the z axis and 
differ by a 90° rotation in a similar way as discussed for the pitchfork bifurcation of 



the excited state in section |3.1[ The resulting branches of the pitchfork bifurcation P2 
in the ground state branch can be seen as functions of the scattering length a — ap2 
and the asymmetry parameter s in figure [2|c) and (d), respectively. Up to a change 
in sign of the splittings AE^f and Aw the results qualitatively agree with those of the 
bifurcation PI of the excited state in figure |2]^a)-(b). The eigenvalues and eigenvectors of 
all three states coincide at the critical scattering length ap2 = —0.0036 of the pitchfork 
bifurcation. 

3.3. Linear model with non-Hermitian matrices 

The bifurcations and exceptional points in dipolar condensates result from the 
nonlinearity of the GPE ([T]). Nevertheless, we now introduce a linear model with non- 
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Hermitian matrices which can reproduce the pitchfork bifurcations shown in figure |2] 
and the level splittings of both the mean-field energies AE^fj and the values Awj as 



defined in equations (19) and (20), respectively. Furthermore, the linear matrix model 



can describe the observed structures when an exceptional point related to one of the 
pitchfork bifurcations is encircled in either the asymmetry parameter s or the scattering 
length a as discussed in section [4], and is even valid for the two parameter perturbations 
discussed in section I4.2[ 

The values of Awj are obtained as eigenvalues of the non-Hermitian matrix 

/ 1 \ 



Q^ 



1 






(26) 



where s is the asymmetry parameter, and 

d = Ca{a - Ocr) (27) 

is the rescaled and shifted scattering length so that the critical value of the bifurcation 



is at a = 0. Equation (20) ensures that the matrix is traceless, and and Ca are real 



parameters, which depend on the values 7 and A of the trap potential and are determined 
by a least-squares fit using the data obtained with the TDVP as described in section [2] 
For a = and s = the matrix has the Jordan form of an EPS with eigenvalue 
A = (see section |4]). The solid lines for the real and imaginary part of Aw in figure |2] 
present the eigenvalues of the matrix with Ca = 26.180 and = —1.0780 x 10""^ for 
the pitchfork bifurcation PI and = 104.46 and c^, = 1.8436 x 10^^ for the pitchfork 
bifurcation P2. For s = one eigenvalue is zero and two eigenvalues are ~ y/a — Ocr, 
and for a = all three eigenvalues are ~ s^^^. Obviously, the results of the matrix 
model perfectly agree with the exact data marked by symbols in figure |2} 

The matrix model for the mean-field energy is somewhat more complicated. For 
the simplicity of the model we require that the matrix elements are low-order power 
functions of d and s, and the eigenvalues are ~ for s = and ~ s^^^ for a = 0. 
Furthermore, we require that for s 7^ both matrices Qy^ and Qe always have the 
same number of degenerate eigenvalues when the two parameters d and s are varied (see 



section 4.2 for more details). The matrix 



Qi 



ce 



d^ 

16as 
9s^ 



9s 2d 



52 



9s 

-2^2 



(2J 



fulfils these requirements. The real coefficients and ce are determined by least- 
squares fits, and read Ca = 26.182 and ce = —3136.3 for the pitchfork bifurcation PI 
and Ca = 101.76 and ce = 4.0582 for the pitchfork bifurcation P2 in figure [2j The 
eigenvalues of the matrix Q e shown as solid lines in the graphs for AE^f agree perfectly 
with the exact data marked by symbols. 
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4. Exceptional points 

The coalescence of two eigenvalues and the corresponding eigenvectors is a feature known 
as "exceptional point" [8 10 , and such an EP2 can appear in systems described by 



non-Hermitian matrices, which depend on at least a two-dimensional parameter space. 
When the exceptional point is encircled in the two-dimensional parameter space the 
eigenvalues show a characteristic feature of a square root branching singularity, which is 
that the two eigenvalues permute after one cycle in the parameter space, and the initial 
configuration of the eigenvalues is obtained only after two cycles in the parameter space. 
The degeneracy of more than two eigenvalues and eigenvectors is possible, in principle, 
however, an EPn in general requires the adjustment of (n^ + n — 2)/2 parameters 11 , 
which implies e.g. that 5 parameters are necessary for an EPS. 

The signatures of coalescing eigenfunctions have been studied by Demange and 
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Graefe 15 for complex non-Hermitian matrices. If these matrices are transformed to 
their Jordan normal form, the type of the exceptional point is given by the size of the 
Jordan block. For example, a block of size three 

(29) 

characterises an EPS. Now a perturbation can be added to Jeps, viz. 

J*P3 = Jep3 + with P e C^^^ . (SO) 

If Q describes a complex path which encircles the exceptional point, e.g. = e"^, 

the behaviour of the eigenvalues depends on properties of the perturbation. If the 
perturbation matrix P is such that the element 

PsiT^O, (31) 



then the typical cubic- root branching singularity expected for an EPS 11 is observed 



i.e. all three eigenvalues and eigenstates permute. However, if the condition (SI) is not 



fulfilled one might observe the permutation between two states, i.e. the signature of a 



square-root branching singularity 15 



In a dipolar BEG with an axisymmetric trap the eigenvalues and the eigenvectors 
of three states coincide in a pitchfork bifurcation by varying only a single parameter, 
viz. the scattering length a of the contact interaction. The coalescence of three states by 
varying only one parameter is certainly related to an underlying high symmetry of the 
system. In this section we investigate in detail the signatures of the exceptional points 
occurring in dipolar condensates. 

4.I. Encircling the exceptional points 

We examine the permutation behaviour of the stationary states when the bifurcation 
points are enclosed by a parameter path. Using the analytic continuation of the GPE 
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-2 X 10~* 2 X 10"* -2 X 10~* 2 x IQ-^ 

Re (Aui) Re {Aw) 



Figure 3. Encircling of the pitchfork bifurcation PI of the excited state (a) in 
the complex scattering length a with the radius = 10^^ and (b) in the complex 
asymmetry parameter s with the radius = 10~^. The permutation of only two 
states in (a) and all three states in (b) is clearly exhibited, (c) and (d): Same results 
for the pitchfork bifurcation P2 of the ground state with Va — 1.5 x 10^^ and = 10^®, 
respectively. The results of the matrix model (solid lines) agree perfectly with the exact 
data (symbols). 



introduced in section 2.1 or the matrix model of section 3.3 the exceptional point can 
be encircled either in the complex continued scattering length a, 

a{(p) = acr + rac"^^ 

or in the complex continued asymmetry parameter s, 

s{ip) = r^e'"^ , 

with the imaginary unit k introduced in section |2.1[ 

The tangent bifurcation marked Tl in figure [l|a) has already been studied 13 



(32) 



(33) 



If the complex scattering length follows the path in equation (32) around the critical 
scattering length oti = —0.019, the two states participating in the bifurcation permute 
in agreement with previous results 13 . This behaviour is typical for an EP2. 

We now investigate the pitchfork bifurcation of the excited state at api = —0.0079. 
When the exceptional point is encircled by varying the scattering length along the path 

the three states behave as visualised in 
y"^. Obviously, the figure 
does not show the permutation of three states, which is typical for a cubic-root branching 
singularity of an EP3. Rather, only the two states emerging in the bifurcation permute, 



of equation (32) with the radius = 10" 
figure ^a) for the expectation values w of the operator 
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which is the scenario discussed in 15 for the special case that the condition (31) is 
not fulfilled. Nevertheless, the typical permutation of three states is observed when the 
asymmetry parameter s is varied along the path of equation (33) with radius = 10~^. 
By following this path the values of the operator w show the permutation plotted in 
figure[3]^b). All three states permute in the way as expected for an EP3 



11 



and discussed 



m 



15 for the generic case that the condition (31) is fulfilled. The different permutation 



behaviour for both control parameters becomes evident in the matrix model introduced 
For s = the matrix element P31 of the perturbation in equation (26) 



m section 



3.3 



I.e. 



vanishes and thus a square root behaviour is expected for two of the eigenvalues 15 
one eigenvalue is constant, Awi = 0, and two eigenvalues follow the paths Aw2,3 
and permute after one circle around the exceptional point. When the exceptional point 
is encircled in the complex s parameter, the condition P31 7^ is fulfilled, and the three 
eigenvalues are ~ e"^*^/^, resulting in the typical permutation of all three states as 
expected for an EP3. 

The pitchfork bifurcation P2 occurring in the ground state for an ansatz with = 6 
coupled Gaussian functions can be analysed in the same way and exhibits a similar 
behaviour as the bifurcation of the excited state. For the path in the scattering length 
(equation (32)) with radius = 1.5 x 10~^ encircling the bifurcation at ap2 = —0.0036 
the values of the operator w show the permutation of two states (see figure [3]^c)). 
However, if the asymmetry parameter s follows the path given in equation (33) with 
all three states permute as shown in figure pi 



Ts = 10- 

As already discussed above, the mean-field energy near the points PI and P2 in 
figure [2] does not show the typical level splittings expected for a pitchfork bifurcation, 
and this is also true when the exceptional points are encircled. The analysis of the matrix 
model Qe for the mean-field energy in equation (28) reveals that two eigenvalues Ai^mf 
follow a path ~ e^'^'^ when the exceptional point is encircled in the scattering length 
a, i.e. one circle in the scattering length results in two loops in the mean-field energy. 
When the exceptional point is encircled in the asymmetry parameter s all three mean- 
field energies follow paths ~ e'*^^/^-''^, which, however, also means a permutation of all 
three states after one circle in s. 



4-2. Two parameter perturbations 

In the previous sections we examined situations where only one of the parameters s or 
a — Oct is nonzero. We now investigate perturbations of the pitchfork bifurcations where 
both parameters are nonzero. Results where either s or a — is set to a constant 
nonzero value and the other parameter is varied are presented in figure |4| In figure |4]^a) 
the symmetry of the trap is broken by setting s = — 10~^. The pitchfork bifurcation 
does no longer exist and only a tangent bifurcation between two states remains, as 
schematically illustrated in figure [l](b). The typical square root behaviour of a tangent 
bifurcation can be clearly seen in the real part of Aw in figure |4f^a), the branch below 
the bifurcation point of the tangent bifurcation belongs to complex eigenvalues and is 
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2 X 10""' 4 X 10~* 

a — ap2 



b) 




-8 X 10"" 8 X 10""' 



Figure 4. Level splittings /S.E and Aw for two parameter perturbations of the 
pitchfork bifurcation P2. In (a) the cylindrical symmetry of the trap is broken with 
s = — 10"'^. In (b) the scattering length is chosen such that a — ap2 = 9.35 x 10"^. 
The symbols mark the results obtained using an ansatz of six Gaussian functions and 
agree perfectly with the solid lines obtained with the linear matrix model. 



obtained only with the analytic continuation of the GPE. In figure |4|^b) the scattering 
length is set to a constant value a = ap2 + 9.35 x 10~^ and the asymmetry parameter s is 
varied. Here, the pitchfork bifurcation is replaced with two tangent bifurcations located 
at s = ±3.6 X 10^^. Three real states exist in the region between the two tangent 
bifurcations, outside that range there are one real and two complex states. 

The symbols in figure |4] mark the results of the GPE obtained with a variational 
approach using six coupled Gaussian functions. They are in excellent agreement with 
the results of the matrix model (26) for the expectation values w = {iplx"^ — y'^\ilj) and 
the matrix model (28) for the mean-field energy. The perfect agreement even for the two 
parameter perturbations is remarkable since both matrix models only use two adjustable 
parameters to describe the level splittings in equations (19) and (20), which have been 
determined solely with the data of the one parameter perturbations. Similar results as 
shown in figure |4] for the two parameter perturbations of the pitchfork bifurcation P2 
of the ground state are also obtained (but not shown) for the pitchfork bifurcation PI 
of the excited state. 

How do the exceptional points related to the pitchfork bifurcations behave under 
the two parameter perturbations? In order to answer this question we follow different 
parameter paths marked CI to C3 in figure |5] either in the complex scattering length 
plane for a constant value s 7^ (see figure [5|a)) or in the complex s plane for a constant 
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Figure 5. Two parameter perturbations of tlie pitclifork bifurcation PI using an 
ansatz with a single Gaussian function and the trap parameters 71 ~ 34000 and Ai = 6. 
In (a) the axial symmetry of the trap is broken with s — 10"'*, and in (b) the scattering 
length is a — api = 2 x 10^^. (c) Illustrations of the composed permutations of the 
states when one, two, or three exceptional points are encircled along the paths CI to 
C3, respectively. 



value a 7^ Ocr (see figure |5[b)). For each path we observe a different permutation 
behaviour of the three states, which can be explained as follows. The colours (or gray 
values) in figure |5]^a) and (b) indicate the minimal distance (absolute value) between the 
mean-field energies of the three states, i.e. a zero value means the coalescence of at least 
two states. In figure |5](a) three points where two states coalesce are revealed. The point 
marked El is a tangent bifurcation point on the real already discussed (see 

figures [l](b) and|4]^a)). The points E2 and E3 are new points in the imaginary half planes 
of the scattering length where different state pairs coalesce. If each point is encircled 
separately they show the permutation behaviour of an EP2 where the two participating 
states permute. Following the path C2 which encircles El and E2 the cyclic permutation 
of all three states is observed. If all points El to E3 are encircled along the path C3, 
we find the same permutation behaviour as for the pitchfork bifurcation of the dipolar 
BEC in an axisymmetric trap, i.e. two states permute. The splitting of the EPS of the 
pitchfork bifurcation into three EP2 for the broken trap symmetry can only be observed 
in the permutation behaviour of the states if a path is chosen which does not include 
all three points El to E3. The points El to E3 in figure [5]^a) merge at a = in the 
limit s — > of the asymmetry parameter. 

The results in the complex s plane for constant scattering length a > Ocr are 
presented in figure |5](b). Here, the pitchfork bifurcation point PI is split into the two 
tangent bifurcation points E4 and E5. If the two points are encircled separately each 
of them shows the permutation of two states typical of an EP2. If they are encircled 
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together along the path C2 in figure [5[b) the cychc permutation of all three states as 
in figure |3](b) is observed, i.e. the typical behaviour of an EPS. As in figure [5|^a) the 
splitting of the EPS into two exceptional points with a square root behaviour can only 
be observed if a path is chosen which does not include both points E4 and E5, and these 
points merge at s = in the limit a — )• a^r- 

The possible permutations of states when one, two, or three exceptional points are 
encircled by a path CI, C2, or CS is illustrated in figure [sj^c). The combinations produce 



the same permutation behaviour for two and three EPs as discussed in 21 



The exceptional points observed in figure [sj^a) and (b) can also be obtained with 
the matrix model in equation (26) or (28). The eigenvalues of the traceless matrices 
Qw or Qe are the roots of the characteristic polynomial xi.^) = + pX + q- Two 
eigenvalues are degenerate when the discriminant D 
the condition 



4p^ + 27g^ vanishes. This yields 



45^ 



27s' 



for the scattering length and the asymmetry parameter. 



(34) 

For a constant value s 



equation (S4) provides one real and two complex solutions for EP2 exceptional points 
in the complex a plane, which agree with the points El to ES marked in figure |5](a). 
For a given value a ^ of the scattering length equation (S4) provides two real or 
two imaginary solutions for s, which coincide with the exceptional points E4 and E5 in 
figure ^h). 



5. Summary 

When dipolar BECs are described within a mean-field theory the GPE exhibits a rich 
variety of nonlinear phenomena including tangent and pitchfork bifurcations of states 
and the occurrence of exceptional points. We have solved the GPE with an extended 
variational approach using coupled Gaussian functions, and presented a method to find 
all states participating in the pitchfork bifurcations, including the complex branches, 
which have been obtained by analytic continuation of the GPE using bicomplex numbers. 
We have analysed in detail the various bifurcations between the states depending on two 
parameters, viz. the scattering length a and the parameter s breaking the axial symmetry 
of the harmonic trap. The origin of the bifurcations is the nonlinearity of the GPE, 
nevertheless, the mean-field energies i?mf and the expectation values w = {tplx"^ — y^\il)) 
of the states participating in the pitchfork bifurcations can be excellently described by 
a linear model with non-Hermitian matrices. 

Both the variational computations and the linear model have been used to 
investigate the properties of the exceptional points. At the bifurcation points of a 
pitchfork bifurcation not only the three eigenvalues but also the eigenvectors coincide, 
indicating the existence of an EPS. However, a different behaviour for the permutation of 
states is observed when the exceptional point is encircled either in the complex continued 
scattering length a or in the asymmetry parameter s. The structures resemble those 
obtained in a linear model using perturbation theory for non-Hermitian operators [15j. 
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For two parameter perturbations the pitchfork bifurcation is spht into either three 
or two tangent bifurcations located in the complex a or s parameter plane, respectively. 
The typical signature of an EP2 is obtained when a single exceptional point is encircled. 
Paths surrounding two or three exceptional points yield the behaviour of combined 
permutations 



21 



The results presented in this article need not be restricted to dipolar condensates. 
Rather, the signatures of exceptional points related to pitchfork bifurcations discussed 
here for dipolar condensates may be generic features of nonlinear systems with pitchfork 
bifurcations, however, further investigations will be necessary to clarify this point. We 
also expect the results of our investigations to arouse the interest of experimentalists 
to search for experimental evidences of the appearance of exceptional points in Bose- 
Einstein condensates. 
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Appendix A. Integrals for the variational ansatz 

The ansatz ([T]) with coupled Gaussians for the wave functions requires the calculation 



of Gaussian integrals to set up the linear system of equations (10). The integrals can 
be computed analytically or with the help of elliptic integrals. With the abbreviations 

^ki = Al + {AlY , af^ = a^' + a'^ with a = x,y,z, (A.l) 

^kl ^ ^fc ^ ^^ly ^ ^kUj ^ ^kl ^ ^ (A^2) 

(A.3) 



klij ij f^i 

klij _ a^_aza^ kiij 

^ ij^klklij ' S 



klij ij f^i 

Ojy az a^ 

ij klij 

Oiy Qiy Ojz 



the integrals read 



7r3/2e-T 



{gV) = V- , (A.4) 



iaW) = 3 : , , (A.5) 



y 



03/2-7'=' 

(aW) = J , , (A.6) 



{gW\9') = . . ! ^ , (A.7) 



3/2„-7'=' 
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{9'\x'V,\g^) 



5/2 



. . , / klij klij klij 

*j=i \l ax ay 'az 



N 



klij 



klij klij 

y az 



yklij 



i 1-1 A I n^^'-i n^^'-i n^^'-i 



klij 



3 / . ' / klij klij klij 

*j=i \l ax ay 'az 



o fciii {'^xl^yRo [l^l, 1^1, l) l) 



{f^yRn {f^l, f^l, l) + 'iu'lKyRx {f^i, fi^, l)) 



(Ka.i?i5 (k^, kJ, l) + 2KxKlRy [kI, kI, l)) ^) 



dux 
Ok, 



4vr5/2 



AT 



klij 



3 . ' / klij klij klij 

*j=i \J ax ay 'az 



2a 



{HxKyRD {kI, kJ, l) - l) 



- ( {fiyRD {kI, 4, l) + 2KlKyRx [nl, Ky, l)) 



- {kxRd {i^l, 4' l) + '^i^xi^lRy {i^l^ 4' l)) 



yklij 



3 . ' / klij klij klij 

*,j=i \ ax 'ay 'az 



2a 



klij 



{KxKyRD {kI, kJ, l) - l) 



dtvx 



(A.8) 
(A.9) 
(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



In equations (A.10)-(A.13) the upper indices at k^'*-^ and k^^^^ have been omitted. The 
elhptic integrals Rd{x, y, z) are given in equation (15), and their derivatives are defined 

as 

Rx{x,y,z) = dxRD{x,y,z) , Ry{x,y,z) = dyRD{x,y,z) . (A.14) 

Integrals of the harmonic trap Vt are easily obtained with the equations given above. 
Integrals not given above are obtained by appropriate permutations of x, y and z. 
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